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Like classical fluids, quantum gases may suffer from hydrodynamic instabilities. Our paper devel- 
ops a quantum version of the classical stability analysis in fluids, the Bogoliubov theory of elementary 
excitations in unstable Bose-Einstein condensates. In unstable condensates the excitation modes 
have complex frequencies. We derive the normalization conditions for unstable modes such that 
they can serve in a mode decomposition of the non-condensed component. Furthermore, we develop 
approximative techniques to determine the spectrum and the mode functions. Finally, we apply our 
theory to a sonic white hole and find that the spectrum of unstable modes is intrinsically discrete. 
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I. INTRODUCTION 

Instabilities may haunt classical as well as quantum flu- 
ids. For example, classical supersonic flows can trigger 
shock waves [[[,§] or moving obstacles in Bose-Einstein 
condensates [|3| can shed vortex pairs Q. In fact, a dy- 
namical instability is at the heart of vortex nucleation in 
rotating condensates j5). In classical fluid mechanics 
the stability of a solution of the equations of motion is 
treated using stability analysis. Assuming a small per- 
turbation of the solution, the equations are linearized in 
the perturbation and the eigenvalues of the linearized 
problem decide the fate of the solution. Complex eigen- 
frequencies with positive imaginary parts indicate insta- 
bilities. In the theory of quantum fluids such as Bose- 
Einstein condensates (3) the equations of motion are lin- 
earized around the mean field to find the elementary ex- 
citations. The ground state of the condensate is, almost 
by definition, stable, yet macroscopic flows of condensed 
atoms may develop instabilities. Here it is important to 
understand how to test for dynamical instabilities and 
how unstable fluctuations evolve. 

In this paper we elaborate a theory of elementary exci- 
tations in unstable Bose-Einstein condensates. Our work 
is primarily inspired by recent proposals ||, 0, |[ |[ 
for generating analogs of black holes using transsonic con- 
densates, but our theoretical concepts and tools may cer- 
tainly find applications in other situations as well [ jLl| . 
Surprisingly little systematic work has been published on 
elementary excitations in unstable condensates, to the 
best of our knowledge, despite the fundamental nature 
of the subject. Inspired by the treatment of instabili- 
ties in quantum fields [ ]l2| , Garay et al. studied a quan- 
tum theory of instabilities in Bose-Einstein condensates 
in an appendix and in a brief book contribution ||. 
Yurovsky |l3}| developed an alternative theory of insta- 
bilities in quantum fluids. Section II of our paper elabo- 
rates on these ideas, starting from the basic concepts of 



elementary excitations in dilute quantum gases H]. Wc 
put forward an economic notation that allows us to derive 
the theory with as little technical effort as possible. Sec- 
tion III addresses two important approximative methods 
to describe unstable excitations analytically. We present 
a brief summary of the frequently applied acoustic ap- 
proximation and develop a version of the WKB approxi- 
mation that can be extended to complex frequencies and 
complex variables. In Section IV we apply all the de- 
veloped concepts and techniques to the analysis of sonic 
horizons, demonstrating so their problem-solving poten- 
tial. 



II. ELEMENTARY EXCITATIONS 

A. Fluctuation field 

Consider a stationary Bose-Einstein condensate of 
atoms with short-range repulsive interactions. Follow- 
ing Fetter jl4| we describe the dynamics of the bosonic 
atom field ip(t, x) using the grand-canonical Hamiltonian 
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Here N abbreviates the total number of particles, 
J ftvbcPx, a conserved quantity, and /i denotes a con- 
stant, the chemical potential. (Because of particle- 
number conservation, both the Hamiltonian H and the 
grand-canonical Hamiltonian fl2.1| ) are equivalent.) We 
assume that most of the atoms constitute a Bose-Einstein 
condensate with macroscopic wave function 
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such that the deviations of if) from the mean field ipo form 
a quantum field <f> of small fluctuations, 



(2.3) 



We expand the grand-canonical Hamiltonian (2.1) up to 
quadratic order in 0, and get 
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with the kinetic term 
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and the condensate flow 
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Here Ho describes the energy of the condensate. The 
Hamiltonian Hi would displace the mean value of the 
fluctuations, when acting on <f>, unless we impose the sta- 
tionary Gross-Pitaevskii equation 
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which minimizes also the Hamiltonian Ho such that 

H = ~ J mc 2 p d 3 x. (2.8) 

Here, and throughout this paper, c denotes the local 
speed of sound, defined by 



mc =gpo- 



(2.9) 



The quadratic Hamiltonian H 2 generates the equations of 
motion of the fluctuation field <j> (Bogoliubov-deGennes 
equations). We found it advantageous to deviate from the 
traditional notation of condensate fluctuations || and to 
combine <f> and cjp in one spinor field 
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Our spinor representation serves as a convenient short- 
hand notation, which does not refer to the spin of the 



atoms of course. In terms of this Bogoliubov spinor (p 
the fluctuation field evolves as 

iHdt<p — Bp , 

B = (T + U — (i + 2mc 2 ) <r z + imc 2 <7 y , 
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where T describes the kinetic energy of the fluctuations. 
Throughout this paper we use the Pauli matrices in their 
standard representation 
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The Bogoliubov-deGennes equation ( 2.11| ) is non- 
Hermitian, because of the anti-Hermitian spinor-mixing 
term imc 2 a y . Therefore, the spectrum of elementary ex- 
citations is not necessarily real. Dynamical in stabi lities 
may emerge. Finally, as a consequence of Eqs. ( [2.11 ), we 
find 

H 2 = jJ (ft(d t 4>) - (dt&)$) d 3 x , (2.13) 

an expression that we need in the mode expansion of the 
Hamiltonian H 2 . 



B. Mode expansion 

As in standard field theories, we expand the Bogoli- 
ubov spinor ip into modes. First we note that tp is invari- 
ant under the conjugation 

= a x tpK (2.14) 

Consequently, tp must have the mode structure 

^ (2.15) 



p = ^ (w v a v + w v a' u ) 



The spinor w v comprises Bogoliubov's familiar u v and v v 
modes j|, 



(2.16) 



and w v denotes the conjugated Bogoliubov spinor, 



W v = a x w u 



such that 
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The mode functions w v are subject to the Bogoliubov- 
deGennes equation 



ihdtw v = Bw v 
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that implies 



ihdfWjj — Bw v . 



(2.20) 



In a field theory, modes are orthonormal with respect to 
an invariant scalar product {w\, 11)2) with 



d t (w 1 ,w 2 ) = 0. 



(2.21) 



in order to guarantee that the d v and a\, are annihilation 
and creation operators. Such a scalar product is 



[wi,w 2 ) = J w\a z w 2 d 3 x, 



(2.22) 



where w' abbreviates w* T . This scalar product is time- 
invariant, because 



a z B. 
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The a„ and a\ are Bose annihilation and creation oper- 
ators if we require that 

{w v , w v i) = J {u* v u v < - v*,v v ) d 3 x = § vv i , 

(w u ,w„>) = / (v v u vl - u v v u i) d 3 x = . (2.24) 



So far, the mode expansion ( 2.1 5| ) is fairly general. Now 
consider single-frequency modes, 



(2.25) 



We obtai n from the invariance (2.21) of the scalar prod- 
uct (|2.22D 



= id t {wi,w 2 ) = (^2 _ Qi)(wi,w 2 ) ■ 
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Consequently, single- frequency modes are only normaliz- 
able when their frequencies f2„ are real 10 v |{L4| . In this 
case we obtain the well-known mode expansion of the 
Hamiltonian H 2 [ pT[ in terms of independent harmonic 
oscillators, 

H 2 = J2 h ^ (*t*> ~ J \v»\ 2 d 3 xj . (2.27) 

Single-frequency Bogoliubov spinors with complex f2, in- 
dicating instabilities, have zero norm. However, this fact 
does of course not prohibit the existence of instabilities. 
It only means that we must not directly employ such 
spinors as modes. 



C. Unstable condensates 

Suppose that the frequency of a Bogoliubov spinor w 
is complex, 



Q = u + ij . 



Because the Bogoliubov-deGennes equation (2.19) is non- 
Hermitian, each eigenfrequency corresponds to a left and 
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a right eigenfunction, here denoted by u>+ and wi, re- 
spectively, 

Bw+ = hnw + , w\B = ffl.w\. (2.29) 
As a consequence of the property ( 2.23| ) we find 

B{a z wi) = a z B^wi = tel*(a z wi) . (2.30) 



Therefore, a z u>i is the Bogoliubov spinor with the com- 
plex conjugated eigenfrequency of w+. Complex frequen- 
cies of elementary excitations occur in conjugated pairs, 
reflectin g th e Hcrmiticity of the grand-canonical Hamil- 
tonian (2.1). The spinor conjugate W+ corresponds triv- 
ially to the frequency — D,*, whereas the frequency — Q is 
associated with the spinor 



w_ = <J z Wl . 



(2.31) 



Consider the scalar product of the modes with ±fi n fre- 
quencies that are labeled by the subscripts ±n, 

= id t (w-n,w+n>) = ("fin + £ln>)(W- n , w +n >) . (2.32) 
Consequently, we can require 

(W-n,W+ n >) = <W (2.33) 

by choosing the appropriate overlap between the left and 
right eigenstates of B. In Section IV we use this orthog- 
onality condition to find the unstable elementary excita- 
tion of a sonic horizon. 

Single-frequency Bogoliubov spinors with complex fl 
must not represent modes per se, yet nothing prevents 
us from combining two or more of such spinors to form 
non-stationary modes. A simple choice is 



satisfying the relations 



(2.34) 



(W± n , W± n >) = 5 nn , , (W± n ,W Tn >) = 

(W±„,W±„0 = 0, (W±„,W T „0 = 0. (2.35) 

Therefore, the W± n are perfectly suitable as Bogoliubov 
modes. We expand the fluctuation field 4> in terms of the 
U±,V± modes, 



i> n = U +n a +n + V+ n a J +n + U- n a- n + Vl n aLp,.36) 



We obtain from Eqs. ([2.25D and (|2.28| ) 



id t U± r , 
id t V± r , 

and, consequently, 



±U! n U± n - ij n V^ n , 
±U n V± n - IJvU^n ■ 



(2.37) 
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id t (f> n = (uj n a +n - i7„at„J U +n - (u n a+ n + i7„a_„J V* n - (uj n a- n + ilnO+n) U- n + (uj n aL n - i"f n a+n) V*„ . 

(2.38) 

We insert this result and the expansion ( 2.36 ) into formula ( 2.13 ), use the orthonormality relations ( 2.35 ), and get 
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Due to the instability of the condensate, pairs of elemen- 
tary excitations are spontaneously generated at the rates 
7„ , and so the non-condensed part grows at the expense 
of the condensate. Of course, the Hamiltonian H2 de- 
scribes the correct dynamics only for short times, as long 
as the growing excitations are still small compared with 
the condensate. Furthermore, the backaction of the non- 
condensed part onto the condensate ought to be taken 
into account, affecting the growth rates 7„ and the fre- 
quencies uj n . The instability causes the condensate to 
dissolve. Nevertheless, the atoms may settle afterwards 
to constitute a new condensate with a stable macroscopic 
wave function, as happens in vortex nucleation |s) . 

III. APPROXIMATIVE METHODS 
A. Acoustic approximation 

Frequently, approximative methods provide the tools 
to find analytic results that capture the essential physics 
of elementary excitations. The best known example is the 
excitation spectrum of a condensate in a harmonic trap 
flBI . Here the excitations of the condensate have been 
calculated in hydrodynamic or, as we would prefer to call 
it, acoustic approximation. (Elementary excitations are 
sound waves within the validity of the approximation.) 
Furthermore, sound waves in moving Bosc-Einstein con- 
densates propagate in the same way as mass-less waves in 
a sufficiently large class of curved space-time structures 
[PI, |7[ p|, 16[. In Section IV we use this connection to 
analyze the instabilities of a sonic horizon. 

Let us briefly summarize the main as pects of the 
acoustic approximation. Given a solution (2.16) of the 
Bogoliubov-deGennes equation ( [2.19] ), the function 

(3.1) 



V> = ip a + e li, °(u + v*) 
solves the time-dependent Gross-Pitaevskii equation 
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as long as u and v are small. We represent tp in terms of 
the particle density p and the phase S, 

4< = ^fpe %S i p = po + p a , S = S + s. (3.3) 



Neglecting the quantum potential ?i 2 (V 2 v 73)/(2my / p) in 
the Gross-Pitaevskii equation ( |3.2| ), we recover the equa- 
tion of continuity and the Bernoulli equation, 
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Assuming that po and 5o satisfy Eq. (|]j) and linearizing 
in p s and s gives 



d t p s = -V- [p s u + po— Vs 

m 



Hpo 



(dt + u- V)s . 



(3.6) 
(3.7) 



where c denotes the speed of sound (2J) and u describes 
the flow (/i/m)V5o. Substituting the expression for p s 
produces the wave equation for sound in irrotational flu- 
ids @,0 



0. 



fl xv = PO ( 1 U 

J c 2 V u -c 2 l + U(g> u 

d v = (d t ,V), 



(3.8) 



the central argument in the analogy between sound in 
moving media and waves in general relativity [l6|| . 
Here we have used a relativistic notation with p and v 
referring to space-time coordinates (not to the chemical 
potential of course). 

To see how p s and s are rel ated to t he B ogoliubov 
spinor we compare tjj 2 of Eqs. (3T) and ( |3.3| ) to linear 
order in p s , s, u and v, and get 



u + v 



(3.9) 



Consequently, a single-frequency sound wave appears as 
the excitation mode 



2mc 2 
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(-in + u • V) + 1 J iae- int , 
(-iil + u • V) - 1 ] iae- mt (3.10) 
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where 



int 



(3.11) 



B. WKB approximation 



Another important approximative method to analyze 
elementary excitations is the WKB approximation Jl8| , 
the equivalent of semiclassical wave mechanics or geomet- 
rical optics (geometrical acoustics). Frequently, the im- 
portant features of waves are determined by the turning 
points of rays on the complex plane or by branch points 
of the momentum. Here we develop a modification of the 
WKB approximation for elementary excitations that 
can be analytically continued to complex variables and 
complex frequencies. 

We assume an effectively one-dimensional model with 
spatial coordinate z and flow u. We express the Bogoli- 
ubov spinor as 



w = (wo + hwi + . . . ) exp 



pdz — Et 



(3.12) 



in terms of the semiclassical mom entum p and the en- 
ergy E. We insert th e ansa tz ( 3.12 ) into the Bogoliubov- 
deGennes equation (2.19) and expand the result into 
powers in h. In zeroth order we obtain 



B w = 
Bo = 



Ew Q , 
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The determinant of Bq — Et vanishes when p satisfies 
the Hamilton-Jacobi equation 
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(E — up) = m c 



with 



E = U H u - (1+ 2mc 2 



(3.14) 



(3.15) 



The components uo and vn of the envelope wq are linearly 
dependent, because wn is an eigenstate of B . We get 



u = vuo, 
1 

V = 
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Before we proceed, let us see how the Hamilton-Jacobi 
equation is related to th e a coustic approximation. The 
hydro dynamical model (3.4) presumes small local vari- 
ations in the density and in the flow such that we can 
describe the condensate locally. Let us assume that the 
chemical potential \x corresponds reasonably accurately 



to the local energy of the condensate, mc 2 + U + 
such that 



En — mc 



(3.17) 



Then, for \p\ <C mc, i.e. for sound with wavelengths 
much larger than the healing length fi/(\/2mc) ||], we 
get the dispersion relation of sound in moving media 



(E -up) 2 = c 2 p 2 



(3.18) 



Additionally, in the relationship (3.16) between the enve- 
lope components we ignore the p /(2m) term within the 
acoustic approximation. We get 



E - 
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(3.19) 



which represents the lin earized Bernoulli equation (3.7) 
expressed in terms (3.9) of the Bogoliubov-spinor com- 
ponents. 

To first order in h we obtain from the Bogoliubov- 
deGennes equation ( 2.19| ) and from the ansatz ( 3.1 2| ) 



= Biw + a z (B Q - Et)wi , 

B 1 = ^—((2pd z +p')t + m(2ud z +u')a z )('i.20) 
2im \ / 

For complex frequencies or complex z values the matrix 
&z(Bn — EH) is symmetric but not necessarily Hermi- 
tian. Since wq is the eigenvector of <j z (Bq — Et) with 
zero eigenvalue, B\Wn must be orthogonal on wq with re- 
spect to the scalar product wfut2- Note that this scalar 
product does not involve complex conjugation. We find 
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= 2^ 9 ^{( u o + v o)P + m i u o - v o) u ) 
which gives the continuity relation 



d z (ul - vl)v = 



with the velocity 



1 + rf p 



1 



rj 2 m 



dH 

dp 



in terms of the semiclassical Hamiltonian 



E = H = up± \j [ +E Q ) -m 2 c 4 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



The continuity relation ( |3.22| ) shows how the spinor am- 
plitudes are connected on the complex plane. For real 
energies and real coordinates (|wo| 2 — |^o| 2 )^ is exactly 
conserved, describing a stationary quasiparticle flux |L8j. 
In this case our result agrees with (|wo| 2 — |uo| 2 )f, apart 
from an arbitrary phase. On the other hand, for com- 
plex energies or complex coordinates (v.q — Vq)v does not 
correspond to an exact conservation law in general, but 
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still remains a constant within the validity of the WKB 
approximation. The advantage of our result is that it can 
be analytically continued such that 



-4n 



(3.25) 



with a constant A Q as long as z does not reach the vicinity 
of a turning point or crosses a Stokes line in the complex 
plane fil. 



C. Turning points 

At a turning point the velocity v vanishes and, conse- 
quently Uq — Vq diverges such that the WKB approxima- 
tion is no longer valid in its vicinity. Turning points are 
the origins of Stokes lines where the WKB solutions are 
discontinuously connected Q. 

Assume that the condensate is nearly uniform around 
a turning point such that the speed of sound is approxi- 
mately constant and that the locality condition ( 3.17 ) is 
satisfied. In this case the turning point does not lie at the 
edge of the condensate, as it is the case for oscillations in 
harmonic traps Jl5[ ] . Let us find out whether and where 
such turning points exist. Our conditions imply that the 
momentum p depends on z only through the flow u{z). 
Let us turn matters around and regard z as a function of 
u and u as a function of p. We get from the expression 
(3.24) of the quasiclassical Hamiltonian 



dH d . . du 

v = — — = u + — {kj — up) = —p-^r- ■ 



dp 



dp 



dp 



(3.26) 



The defini tion ( 3.16 ) of r\ and t he H amilton- Jacobi equa- 
tion (3.14) with the condition (3.17) implies the relation 



because the energy of the elementary excitation oug ht to 
be much smaller than the condensate's energy (|2.8|). We 
expand 770 in powers of e 1 / 3 , 



Vo = xo + xie 1/3 + x 2 e* /L ' + x%e 



2/3 



r 3/3 



and find the coefficients 



Xq = -1 , Xi = V-4, X2= , X3 

Xi 



(3.31) 



(3.32) 



Close to a turning point the flow profile depends quadrat- 
ically on the momentum 



U — Uq 



d 2 u 
dp 2 



ip - Pa) 2 



(3.33) 



This relation shows that a turning point is a branch point 
for the semiclassical momentum p and it also specifies 
the onset of Stokes lines [19|, defined as the lines where 
the differences between the phases J pdz of the two p 
branches is purely imaginary. Here one of the waves is 
exponentially small compared with the other. Crossing 
a Stokes line connects waves in a discontinuous yet pre- 
cisely defined w ay [Jl9| . We obtain the coeff icients from 
Eqs. ( 3.27 ) and (|3.30| ), using the series ( 3.31 ), 
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(3.34) 



Finally, w e det ermine the velocity uq at a turning point 
from Eq. (|||). We find 

/ - 2/3 

(3.35) 



V 



(3.27) 



We substitute this result for p 2 in Eq. ( 3.16 ), solve for 
up, and get 



1 



2,/ 



with 



E 



mc 



2 ' 



(3.28) 



(3.29) 



We differentiate Eq. ( 3.28|) w ith respect to the momen- 
tum p, utilize relation Q3.27 ) and its momentum deriva- 
tive, and arrive at 



du _ 1 + (3 - 2e)r) + (3 + 2e)r] 2 
dp ~ 2m(r] - 1)(t? + l) 2 



(3.30) 



The turning points correspond to the zeros 770 of the nu- 
merator. We note that e is naturally a small number, 



Here ty—I refers to the three cubic roots of —1, gener- 
ating three turning points in the complex plane when u 
approaches c. In a mostly uniform condensate the turn- 
ing points of elementary excitation are close to transsonic 
regions where the condensate flow transcends the speed 
of sound. Such a region forms a sonic horizon. 



IV. SONIC HORIZONS 
A. Model 

Consider a stationary transsonic medium i.e. a 
medium with a spatially non-uniform flow that varies 
from subsonic to supersonic speed. One would expect 
that beyond the interface where the flow exceeds c sound 
waves are swept away such that no sound from the su- 
personic zone can return to the subsonic region. This 
transsonic interface serves as the sonic equivalent of a 
black-hole horizon. || [t|, |8|, On the other hand, 

the interface where the flow settles from supersonic to 
subsonic speed forms the horizon of a sonic white hole 
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H 0, ||, an object that no sound wave can enter from 
outside. Close to the horizons, sound waves propagat- 
ing against the current freeze, and their wavelengths are 
dramatically reduced. 

Transsonic Bosc-Einstein condensates offer great 
prospects || @, ||, |l0) for demonstrating the quantum 
effects of event horizons [^0[ ^f| . Such effects and the sta- 
bility of the condensate depend on the behaviour of the 
condensate close to the horizons. Let us thus focus on 
the physics in the vicinity of a black- or white-hole hori- 
zon. In this case we can use the simple one-dimensional 
model 

u = —c + az . (4-1) 

Here z denotes the spatial coordinate orthogonal to the 
horizon at z = 0, a characterizes the surface gravity or, in 
our acoustic analog, the gradient of the transsonic flow, 
and po and c are assumed to be constant. Strictly speak- 
ing, we should complement the flow component (4.1) in 
the z direction by appropriate components in the x and 
y directions, in order to obey the continuity of the flow. 
But as long as we focus on effects on length scales smaller 
than \c/a\ we can ignore the other dimensions of the fluid. 
Depending on the sign of a, two cases emerge. When 
the velocity gradient is positive we are considering the 
horizon of a sonic black hole. When a is negative the 
horizon refers to a white hole || ^[ ||| . In a typical alkali 
Bose-Einstein condensate without exploitation of Fesh- 
bach resonances the speed of sound c is in the order of 
1 mm/s. The transsonic velocity gradient should be small 
compared with the healing length || 



M£«c, £ 



nic 



(4.2) 



which guarantees that the Hawking energy k B T H is much 
smaller than the energy of the condensate, 



h\a\ n 
Ztt 



(4.3) 




FIG. 1: Analytic structure of a sonic horizon in a Bose- 
Einstein condensate. The figure shows the three turning 
points in the complex z plane around the horizon at the 
origin. The points are connected by two branch cuts where 
the wavenumber k of acoustic elementary excitations is ele- 
vated to other solutions of the dispersion relation. The circle 
around the turning points roughly indicates the place where 
k approaches the asymptotics Q/(az) that is characteristic of 
waves at a horizon |2l| . The third branch cut, connecting one 
of the turning points to infinity, is the cut between the two 
trans-acoustic branches with the asymptotics (4.11). 



We use the relationship ( 3.10 ) to find the Bogoliubov 
spinor within the acoustic approximation, 



w 



A 



where 



/ n 

V 2az 



A 



h 

mc 



-Sq- 



(4.7) 



(4.8) 



Being in the vicinity of the horizon and having the linear 
velocity profile (U.fl) presumes that 



1*1 « 



(4.4) 



The energies of elementary excitations should be suffi- 
ciently smaller than mc 2 , which implies that their di- 
mensionless energy parameter s defined in Eq. ( |3.29 ) is 
small. The excitations are sound waves for low wavenum- 
bers and for z far away from the turning points, 



Given the linear velocity profile fl4.f| ), we solve the wave 
equation (^3) exactly, and get 



s = s ( iz ±in / a e - mt - iz^ in '' a e in 



■')■ 



(4.6) 



The approximation (4.7) is restricted to complex z vari- 
ables outside the trans-acoustic region indicated in Fig. 
I. Close to the trans-acoustic zone the wavelength of 
elementary excitations is dramatically reduced. In this 
regime we can use the WKB approximation. In the re- 
gion where both the acoustic and the WKB approxima- 
tion are applicable we represent the Bogoliubov spinor as 



Uo 

v 

with the wavenumber 



exp 



kdz — iflt 



az 



(4.9) 



(4.10) 



If the acoustic approximation were universally valid the 
wavenumber k would approach a singularity at the hori- 
zon where, consequently, the wavelength of sound would 



8 




FIG. 2: Wavenumbers k of elementary excitations around a sonic white-hole horizon, analytically continued on the complex 
plane. The figure shows three roots of the dispersion relation [ft 2 fc 2 /(2m) + mc 2 ] 2 — h 2 [Q — (— c + az)k] 2 — m 2 c 4 for = 
O.li (mc 2 /fi) and a — —0.5 (mc 2 /ft), illustrating the branch cuts of k. The top row displays the wavenumber of a sound wave 
that propagates against the current. The picture indicates the characteristic Q/(az) asymptotics away from the branch points. 
The two lower rows display two trans-acoustic branches of k. The fourth root of the dispersion relation is not shown, because 
it corresponds to the trivial case of sound waves that propagate with the flow. 



shrink to zero. Apparently, Nature tends to prevent 
such extreme behaviour p^ |. In our case, the singular- 
ity is split into three branch points, the turning points. 
The branch cuts between the points connect the acous- 
tic branch to two other trans-acoustic branches with 
wavenumbers so high that the branches do not repre- 



sent sound waves subject to the wave eq uatio n (p.q). 
We solve the Hamilton-Jacobi equation (3.14) in the 
|f2| <C \uk\ limit and get the asymptotics for the trans- 
acoustic branches E3l 



k 



±2 — y/u 2 /c 2 



1 + 



flu 



(4.11) 
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Figure 2 illustrates the momenta k of the three branches 
in the complex z plane. Finally, we note that the fourth 
solution of the Hamilton- Jacobi equation (3.14) corre- 
sponds to sound waves that are swept away by the cur- 
rent, a less interesting case. 



B. Stokes phenomenon 

It is tempting to assume that we could employ the 
acoustic asymptotics Q/(az) of the wavenumber k on 
the entire complex plane, as long as z is sufficiently far 
away from the turning points at the horizon. In this 
case, however, the phase J kdz becomes logarithmic and 
hence multivalued, whereas the true Bogoliubov mode 
function is singlevalued. The Stokes phenomenon ]1E(] 
resolves this conflict by connecting the acoustic modes 
to trans-acoustic ones on the upper or the lower half 
plane, giving rise to connection formulas p3[ ] that de- 
scribe mode conversion p4| . Each turning point zq con- 
nects two branches of WKB wavenumbers, as we see from 



a.b / 



Eq. (3.33), and each z is origin of three Stokes lines. At 
a Stokes line the phase difference between the two con- 
nected branches is purely imaginary. Consequently, one 
of the WKB waves exponentially exceeds the other such 
that the smaller cannot be resolved within the WKB ap- 
proximation, if the larger wave is present. In general, the 
Bogoliubov spinor is a superposition of the four funda- 
mental solutions that correspond to the four branches of 
the WKB wavenumbers. 



w 



c A w. + c R w R + c n w n + c n w r 



(4.12) 



When crossing a Stokes line, the exponentially sup- 
pressed solution may gain an additional component that 
is proportional to the coefficient of the exponentially en- 
hanced solution. If we wish to construct a Bogoliubov 
spinor where only the exponentially smaller component 
exits in the vicinity of a Stokes line we must put the co- 
efficient of the larger one to zero. Figure 3 shows the 
Stokes lines of Boboliubov modes at a sonic white hole 
with purely imaginary frequency, a case of importance in 
Sec. IV C. The pairs of letters indicate which branches are 
connected by the lines, and the first letter identifies the 
exponentially dominant branch. The picture shows that 
with the choice of branch cuts made we can construct a 
Bogoliubov mode that is acoustic on the lower half plane, 
by demanding that c c vanishes at the C,A Stokes line of 
the lowest turning point and that c B is zero on the lower 
half plane. Since c D is not connected to the other three 
branches we can put the coefficient to zero throughout 
the complex plane. Such a strategy is not possible for a 
sonic black hole, see Fig. 4, because here the order of ex- 
ponentially dominant or suppressed waves is exchanged. 
Moreover, we cannot find a solution that is purely acous- 
tic on the upper half plane, because of the branch cut 
between the central and the highest turning point in Fig. 
4. We can of course alter the arrangements of branch 
cuts, such that the cut between the two trans-acoustic 




A,C 



A,B 



C,A 



C,B 



FIG. 3: Stokes lines of elementary excitations at a sonic white 
hole with purely imaginary frequencies (dotted lines). The 
pairs of letters indicate which branches of the superposition 
(4.11) are connected by the lines. The first letter of each pair 
identifies the exponentially dominant branch. 



C,A 




C,A 



B,C 



A,C 



B,A 



FIG. 4: Stokes lines of elementary excitations at a sonic black 
hole with purely imaginary frequencies analogous to Fig. 3. 



branches B and C points in other directions. Choosing 
different branch cuts does determine in which half plane 
the excitation wave may be acoustic, but it does not influ- 
ence whether such an acoustic behavior is possible at all 
on either the upper or the lower half plane. Therefore, at 
a sonic white hole unstable elementary excitations with 
the acoustic asymptotics (4.7) exists on one of the com- 



plex half planes, whereas the unstable modes of black 
holes, if any, are always trans-acoustic. 
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C. Instabilities 

Consider the unstable elementary excitations of a sonic 
white hole, assuming the acoustic asymptotics fl4.7| ) with 
complex frequencies ui + ij to b e va lid on one of the com- 
plex half planes. Excitations (4.7) with negative 7 are 
localized near the horizon and are attenuated in time, 
whereas acoustic excitations with positive 7 would grow 
in space and time. Therefore, the truly unstable modes 
are trans-acoustic. Nevertheless, we can use the proper- 
ties of the attenuated modes to determine the spectrum 
of unstable modes, because of the spectral symmetry of 
the Bogoliubov-deGennes equations proven in Sec. II C. 

Bogoliubov excitations with complex frequencies have 
zero norm and should satisfy the orthogonality relation 
(2.33). First we calculate the scalar product (101,103) 
at time r = 0, to find a condition for zero norm. We 
assume that the scalar product ( 2.22| ) is dominated by 
the acoustic region where \zq\ <C \z\ <C \c/a\, 



(wi,w 2 ) = J w{a z w 2 dz 

- \A\ 2 
= \A\ 2 



, mc 
, mc 



2w 



h 

x / — exp 

& z 



f e ±2nw/a + A 



.Ul UJ2 J dz , 413 x 

a 



Here we have connected the two acoustic regions on the 
upper or lower half plane, respectively, indicated by the 
± sign, circumventing the trans-acoustic region close to 
the horizon at z = 0. The scalar product vanishes if 



UJ = 0. 



(4.14) 



Therefore, within the validity range of our approxima- 
tions, the unstable elementary excitations have purely 
imaginary frequencies 7. Let us examine the orthogonal- 
ity condition (p.33|). We represent 7 as 



7 = 2na . 



(4.15) 



and we deform the integration contour such that it cir- 
cumvents the trans-acoustic zone on a large semi-circle 
with radium r from below. 



(w- n ,W+n') = j {v_ n u +n , - U- n V +n >) dz 

~ A- n A™ f 2 -^±^z^dz, 
h J iz 



L ' 



x / e 2i{n - n ' )e d6 . 



(4.16) 



The scalar product (w- n , w +n i) vanishes for n 7^ n! if 
the n are integers. The spectrum of unstable elementary 



excitations consists of discrete and equally spaced points 
on the imaginary frequency axis. 

An interesting proposal for a sonic hole [|[ @, |§| involves 
a toroidal condensate that flows through a constriction 
where it exceeds the speed of sound and that then, after 
the constriction, settles to subsonic speed. Garay et al. 
H 0, H found that the condensate is unstable only in 
narrow "instability fingers" in the parameter space used. 
Our analysis indicates that the instabilities are generated 
when the excitations of the toroidal condensate match the 
discrete imaginary resonances of the white hole. This 
would explain the narrowness of the instability fingers 
H 0, H). Such instabilities are enhanced by the "las- 
ing" effect of the black-hole white-hole pair J24|. The 
white hole generates a hydrodynamic instability that is 
resonantly enhanced by the pair of horizons where ele- 
mentary excitations can bounce back and fourth. Sonic 
black holes can be stabilized by employing the equiva- 
lent of a Laval nozzle || that converts a subsonic flow 
to a supersonic one without causing turbulence (as in a 
rocket engine) . Our theory indicates that white holes are 
intrinsically unstable ^5|, generating breakdown shocks 

I- 



SUMMARY 



Unstable Bose-Einstein condensates develop elemen- 
tary excitations with complex frequencies. Such excita- 
tions have zero norm and are subject to orthogonality 
relations between pairs of excitations with opposite fre- 
quencies. We elaborated the general theory of unstable 
elementary excitations and of two important approxi- 
mate methods to analyze their behaviour, the acoustic 
and the WKB approximations. Applying these tech- 
niques, we showed that sonic white holes in Bose-Einstein 
condensates give rise to a discrete spectrum of instabili- 
ties, which may explain the remarkable stability of sonic 
holes in toroidal traps |^, @, § ■ 
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